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Abstract. The magnetic flux freezing theorem is a basic principle of ideal 
magnetohydrodynamics (MHD), a commonly used approximation to describe the 
aspects of astrophysical and laboratory plasmas. The theorem states that the magnetic 
flux — the integral of magnetic field penetrating a surface — is conserved in time as that 
surface is distorted in in time by fluid motions. Pedagogues of MHD commonly derive 
flux freezing without showing how to take the material derivative of a general flux 
integral and/or assuming a vanishing field divergence from the outset. Here I avoid 
these shortcomings and derive flux freezing by direct differentiation, explicitly using a 
Jacobian to transform between the evolving field-penetrating surface at different times. 
The approach is instructive for its generality and helps elucidate the role of magnetic 
monopoles in breaking flux freezing. The paucity of appearances of this derivation in 
standard MHD texts suggests that its pedagogic value is underappreciated. 
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1. Introduction 

Magnetohydrodynamics (MHD) is the simplest generalization of hydrodynamics for a 
sufficiently ionized collisional plasma [TJ[2j[3]. The relative motion between positive and 
negative charge carriers creates currents which can sustain magnetic fields and electric 
fields are generated by induction. Charge separation and plasma oscillations are assumed 
to occur on small enough spatial and temporal scales that the plasma is considered to be 
neutral on macroscopic scales of interest. As in hydrodynamics, a high rate of particle 
interactions ensures that deviations from Maxwellian velocity distributions are small. 

Non-relativistic MHD limit has been a mainstay of theoretical astrophysics since 
most of the material inside stars and between them is composed of non-relativistic 
magnetized plasma and is commonly treated in the MHD approximation. The limit is 
also widely used in approximating the dynamics and stability of fusion device plasmas. 
Students in astrophysics and plasma physics are typically exposed to MHD either in 
advanced undergraduate or graduate courses. 

The solution of physical problems in MHD requires equations for mass conservation 
equation, momentum conservation, energy evolution equation and the magnetic 
induction equation. The latter is the subject of the present paper and is given in 



where B is the magnetic field, V is the plasma velocity, and vm = ¥~ is the magnetic 
diffusivity in terms of the resistivity rj and speed of light c. The form of Eqn (1) is 
identical to that of vorticity evolution in incompressible gravitational hydrodynamics if 
B is replaced by vorticity weVxV and um is replaced by the kinematic viscosity v. 

Eq ([1]) is derived by starting with Faraday's law ^ = cVxE, where E is the 
electric field. Eliminating E in terms of B and J is then accomplished by use of Ohm's 
law, E + V x B = 77 J, where J is the current density. Use of the non-relativistic 
Amperes law VxB = then relates J and B. The aforementioned Ohm's law 
is derived by subtracting the separate momentum density equations for positive and 
negative charge carriers. The resistive term arises in Ohm's law when there are finite 
but small deviations from Maxwellian distributions of the charge carriers. That these 
deviations are assumed to be weak highlights that the relevance of MHD when many 
collisions between charged particles occurs over dynamical times of interest. Collisionless 
plasmas have more complicated Ohm's laws. 

The induction equation of ideal MHD corresponds to Eq. (CQ) when the Vm term 
is neglected. The ratio of magnitudes of the second term to the third term in Eq. (pQ) 
can be approximated by the magnetic Reynolds number Rm = ^rS where V is the 
velocity magnitude and L is the characteristic gradient scale of velocity or magnetic 
field. When Rm » 1, the resistive term is often ignored and ideal MHD assumed. 
Actually, many astrophysical plasmas are turbulent which implies a spectrum of eddies 
of different scales and energies. Such flows always have a microphysical dissipation 



CGS units by [3] 
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scale at which Rm = 1 so in practice one must think about Rm as a scale dependent 
quantity. Subtleties involved with applying the limit of ideal MHD to turbulent flows 
are discussed in jl] and [5]. For the laminar (non-turbulent) case, the limit Rm >> 1 
leads more straightforwardly to the ideal MHD approximation and that is the focus of 
this paper. 

The ideal MHD limit is indeed the limit for which "magnetic flux freezing" or 
Alfven's theorem holds. This theorem (analogous to the Kelvin circulation theorem of 
ideal hydrodynamics [El [7] in which vorticity flux is frozen) states that the magnetic 
flux through a material surface is conserved even as velocity flows distort that surface 
in time. Mathematically this means 



Dt 



0, (2) 



where D/Dt indicates the material or Lagrangian derivative + V(x, t) ■ V for space 
and time dependent flow velocity V(x, t) and magnetic flux = |B ■ dS, where 
the integral represents a surface integral over an open surface. The rest of this paper 
addresses derivations of Eq. (T5j) . 

In section 2, I derive Eq. (T5]) by a formal material derivative of the flux integral. In 
section 3, I compare this derivation to other derivations commonly found in standard 
texts. In section 4, I address why the derivation of section 2 facilitates a better physical 
understanding of the role magnetic monopoles play in violating flux freezing compared 
to the other derivations discussed in section 3 . I conclude in section 4. 



2. Material derivative of the flux and derivation of flux freezing 

Consider an open surface within a plasma at t = through which magnetic field lines 
penetrate. The surface has a differential area dSo = daida 2 , where 0\ and a 2 define 
local Cartesian coordinates. Via the distorting action of a smoothly varying time and 
space dependent velocity flow V, this surface evolves to a new surface with differential 
dS at time t > 0. To compute D^B/Dt we must then account for the fact that both B 
and dS can depend on space and time. 

The time-evolved surface measure dS can be related to the initial surface measure by 
a Jacobian transformation. With this in mind, let us define a general flux $q = / Q • dS 
for an arbitrary vector Q, not necessarily divergence free. Then 

where dS is the differential surface area vector at an arbitrary t > and dSo = ||c?So|| 
is the scalar differential surface area of the surface at t — 0. Here J7" is the Jacobian 
vector relating dS to a coordinate transformation of dSo, and has components 

sj- 0x r 0x s , . 

Jq = eqrs do T 1 do T 2 ' [) 
where e qrs is the Levi-Civita symbol and xi, x 2 , x 3 are the local Cartesian coordinates of 
the evolving surface element dS. Three coordinates are required for the evolving surface 
as the surface normal can evolve away from its initial direction. 
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Because the right side of ([3J) now involves an integral over a fixed surface, we can 
take the ■&? inside the integral to obtain 

DJ 



^/Q-^o = /(^Q + Q-^)^o. 



(5) 



We now need an expression for DJ'jDt. Changing the indices q, r, s in (j4]) to k, i,j and 
taking the material derivative gives 

DJ k 



Using 



Dxj_ 
Dt 



Dt 

rr . Dxj 
V J' Dt 

DJk 
Dt 



d(Dxi/Dt) dxj + dxi d(Dxj/Dt) 



V- ^ 



d(T\ da 2 do\ d(72 

... dV^dxrn. nnr | W± - 9V J dx m 
dx m aa\ ' oa 2 ox m d<T 2 



in Eq. (jSJ) then gives 



dVi dx m dxj dxt dVj dx r 



dx m d<j\ da2 d<j\ dx m da 2 



Since e^- = — e^, we can interchange indices i and j in © to obtain 



DJk 
Dt 



eki i dxL 



dx m dXj 



dXj dx r , 



dcri d(J2 da\ d(J2 



kt 3 dx m m 3 qU 1 



(6) 



(7) 



(8) 



where the second equality follows from using 03} multiplied by e kmn and the third 
equality follows from using e kij e mjq = 5 qk 5 mi - 5 qi 5 



Using OHD, along with the definition D/Dt and the fact that (V • VQ)» = V • VQi 
in Cartesian coordinates, we now obtain for (j5J) 



D 
Dt 



Q • JdS 



'dQt 
dt 



+ V • VQi + QiV ■ V - Q • W* JidS . (9) 



Using JidS$ = dSi and the vector identity 

Vx(VxQ) = Q VV-V VQ-QV-V + VV-Q, 
Eq. (jUJ) becomes 
D 



Dt 



jQ-dS = J |^-Vx(VxQ)+VV-q| -dS. 



(10) 



Thus if ^ = Vx (V x Q) and V ■ Q = 0, then -g / Q • dS = 0. Replacing Q with B 
we have the proof of flux freezing by direct differentiation. 

The derivation above is just the surface integral analogue of that used to prove 
Reynolds transport theorem. The latter describes the evolution of a scalar volume 
integral over a material volume that evolves in time from a velocity flow [HI [9l [TO] . 



3. Comparison to other derivations 

Like the derivation above, the approaches in [TTj and [12] do proceed by direct 
differentiation of the flux integral and separately compute DB/dt and DdS/dt by 
considering the infinitesimal evolution of these quantities. However these approaches 
make no explicit mention of the Jacobian so the connection to the basic method of 
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surface integral transformations is not made explicit. They also assume V • B = and 
do not provide a physical interpretation of why V-B^O can violate flux freezing. 

The derivation of the Kelvin circulation theorem in Ref. [7] also proceeds by direct 
differentiation, but uses Stokes' theorem first so that the derivative is taken on the line 
integral of velocity. The fact that vorticity is the curl of the velocity automatically 
ensures that vorticity is divergence-free. Also, if the same approach were applied to 
magnetic flux then the mathematical analogue to the velocity in the proof is the vector 
potential which itself is not a gauge invariant quantity. The derivation of section 1 
avoids use of vector potential for the case that Q = B and carries the divergence term 
to the very end. 

Most noteworthy is that in addition to starting with V ■ B = 0, common MHD 
presentations [H El [3j [THl Ej do not cleanly show how to calculate the material derivative 
of the flux integral. Instead these approaches are characterized by the following: A 
bounded open surface C in the plasma is considered to evolve in a small time St to 
a new surface C by the action of differentiable velocity flows. Because V ■ B = 0, 
Gauss' theorem tells us that the total integrated flux at time t + St through the closed 
surface formed by C, C and the quasi-cylindrical "side" connecting C and C is zero. 
Mathematically, this means 

[v-BdV = = - / B(r,t+St)-dS+ [ B(r, t+St)-dS+ [ B(r,t+5t)-dS.(12) 

J JC JC J side 

The last term is 

/ B(r, t+St)-dS = ( B(r, t+St)-(dlxV St) = f <ft(VxB(r, t+St))-dl, (13) 

J side JC JC 

where d\ is the line element around the boundary of surface C. The differential change 
in magnetic flux through C as it evoles from C to C is 

S$ B = I B(r,t + 6t)-dS- I B(T,t)-dS. (14) 

JC JC 

Using (JTSj) and (fT4|) to replace the last and penultimate terms of (fT2|) respectively, gives 
S<$> t f)B r t c)W t 

^ = / c ar dS+ / c ( VxB ( r - t »" fl =/ c iir dS+ / c Vx ( VxB ( r - i )»' dS -( 15 » 

where the last equality follows from use of Green's theorem and the limit of small St 
has allowed replacement of B( ^^ t+St )- B< ( r < t ) \yy dB_ Finally, by writing SQ/St = D&/Dt 
and using Eq. (CD) for um = 0, Eq. (j2J) obtains. 

By comparison to the method of the previous section, derivations along the lines 
of those which follow Eq. f|T2l) - f|T5l) do not as lucidly separate of the material time 
derivative of the integrand from that of the measure. In addition, since the starting 
point assumes V • B = 0, the reader is not provided with the opportunity to understand 
why V-B^0 can violate flux freezing. 

4. Seeing the role of V ■ B 

The need for clarity on why V ■ Q ^ can violate flux freezing is further evidenced 
by the ambiguity of the derivation of the Kelvin vorticity theorem of Ref. [15J. 
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There (in contrast to the approach of [7]) it is implied that any vector Q satisfying 
^-Vx(VxQ) = obeys flux freezing and that the property that Q is the curl of 
some other function (i.e. that V • Q = 0) is unnecessary to prove flux conservation. But 
this is incorrect, and seemingly results in Ref. [15] from an ambiguity in distinguishing 
the partial and material derivative (compare Eq. 44 of Ref. [15] to Eqn. |5] above). 
As seen in Eq. ffTTj) above, a finite divergence term would violate flux freezing even if 
the former condition is satisfied. Note that the last term in Eq. (TIT]) containing the 
divergence actually cancels the hidden divergence term within the penultimate term 
since — Vx(V x Q) includes a term — VV • Q when expanded with vector identities. 
However, the ideal MHD magnetic induction equation does not involve V • B = in its 
derivation. Thus for Q = B, the first two terms on the right of ( TlT]) cancel, leaving the 
divergence term whose physical meaning I now discuss. 

The divergence term of Eq. ( flTT) represents the net advection of field line divergence 
through the evolving surface. For Q = B, a finite V • B = 47rp m would imply the 
existence of magnetic monopoles of magnetic charge density p m by analogy to Gauss' law 
for electric charge. The magnetic field lines emanating from a magnetic monopole have a 
net magnetic flux through any spherical surface surrounding the monopole. To see that a 
net advection of magnetic monopoles through a surface would change the magnetic flux 
through that surface, first consider the contribution form a single monopole of positive 
magnetic charge which has all field lines directed radially outward. As the monopole 
approaches the surface from one side and passes through to the other, the sign of its 
contribution to the flux through that surface changes. An advection of a net density of 
monopoles of one sign through the surface would then by extension also change the flux 
through the surface with time. 

Note that a V • B term in the flux evolution equation need not be the only 
consequence to MHD in a hypothetical plasma of arbitrarily large magnetic monopole 
densities. In the same way that we derive the standard Ohm's law for MHD by 
subtracting electron and ion momentum density equations, we would also have to derive 
a magnetic Ohm's law by subtracting positive and negative magnetic monopole charge 
density equations. The two Ohm's laws would be coupled. Further study of magnetic 
monopoles is beyond the scope of the current paper. 

5. Conclusion 

Commonly used derivations of magnetic flux freezing tiptoe around showing how to 
take the material derivative of a general flux integral. In addition, derivations which 
begin with V • B = from the outset do not provide the opportunity for a physical 
understanding of why V-B ^ could violate flux freezing. The direct differentiation 
method of section 1 overcomes both of these shortcomings. 
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